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ABSTRACT 

We present a detailed study of the environments of a sample of nine low-power (Fa- 
naroff & Riley type I) radio galaxies, based on new and archival XMM-Newton observations. 
We report new detections of group-scale environments around three radio galaxies, 3C 296, 
NGC 1044 and 3C76.1. As with previous studies, we find that FR-I radio galaxies inhabit 
group environments ranging over nearly two orders of magnitude in bolometric X-ray lumi- 
nosity; however, we find no evidence for a tight relationship between large-scale X-ray envi- 
ronment and radio-source properties such as size, radio luminosity, and axial ratio. This leads 
us to conclude that radio-source evolution cannot be determined entirely by the global prop- 
erties of the hot gas. We confirm earlier work showing that equipartition internal pressures are 
typically lower than the external pressures acting on the radio lobes, so that additional non- 
radiating particles must be present or the lobes must be magnetically dominated. We present 
the first direct observational evidence that entrainment may provide this missing pressure, in 
the form of a relationship between radio-source structure and apparent pressure imbalance. 
Finally, our study provides further support for the presence of an apparent temperature excess 
in radio-loud groups compared to the group population as a whole. For five of eight tem- 
perature excesses, the energy required to inflate the radio lobes is comparable to the energy 
required to produce this excess by heating of the group gas; however, in three cases the current 
radio source appears too weak to produce the temperature excess. It remains unclear whether 
the temperature excess can be directly associated with the current phase of AGN activity, or 
whether it is instead either a signature of previous AGN activity or simply an indicator of the 
particular set of group properties most conducive to the growth of an FR-I radio galaxy. 

Key words: galaxies: active - X-rays: galaxies: clusters 



1 INTRODUCTION 

The dynamics of low-power (FR-I: Fanaro ff"& Rileviri974h . twin- 
jet radio galaxies are thought to be controlled by interactions 
between the initially highly relativistic jets and the surrounding 
hot-gas medium, which act to dec elerate the jets significantly on 
scales of typically a few kpc (e.g. iBic knell 1994; Laing & Bridle 
l2002al) . On scales of hundreds of kpc, low-power radio galax- 
ies exhibit a wide range of morphologies, from narrow ly colli- 
mated plumes to large, rounded lobes (e.g. lLaindll993l) . In or- 
der to understand how such structures are produced, and to in- 
vestigate the evolution of FR-I sources on large scales, it is nec- 
essary to have good constraints on the hot-gas environments that 
confine their radio lobes. ROSAT observations first established that 
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the brightest and most well-studied F R-I radio galaxies usually re- 
side in group-scale env i ronments (e.g. Worrall & Birkinshaw! 19941 : 
iKomossa & Bohringerl [19991: I Worrall & Birkinshaw! l2000h . With 
Chandra and XMM-Newton it is possible to study radio-galaxy en- 
vironments and jet / environment interactions in more deta il (e.g. 
iBlanton et al1l200ll ; ICroston et al]|2003l : iFabian et al]|2003l) ; how- 
ever, while there have been many studies of the apparent impact of 
FR-I sources in clusters, to date there have been few detailed Chan- 
dra and XMM-Newton studies of the environmental properties of 
nearby, well-studied FR-I radio galaxies. The role of environmental 
properties such as group/cluster richness, gas density and tempera- 
ture distribution in determining large-scale radio structure therefore 
remains uncertain. 

Another key uncertainty both in modelling radio-galaxy dy- 
namics and in establishing their energetic impact is the lack of 
constraints on particle content in FR-I radio jets and lobes. In the 
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case of powerful (FR-II) radio galaxies, two lines of evidence sug- 
gest that synchrotron minimum-energy estimates are roughly valid 
in the lobes: (1) for cases where external pressure measurements 
are available, they are in reasonable agreement with the internal 
minimum pressure (e.g. lHardcastle et alj2002l ; ICroston et al.l2004l ; 
iBelsole et alj|2004h ; (2) measurements of inverse-Compton emis- 
sion from the lobes of FR-II radio galaxies suggest that their elec- 
tron energy densities are close to the value for equipartition with 
the magnetic field in t he absence of an energetically important 
proto n population (e.g. ICroston et all l2005bl ; iKataoka & Stawarj 
2005). However, in low-power (FR-I) radio galaxies, it has been 
known for several decades that a relativistic electron population in 
equipartition with the magnetic field cannot provide sufficient pres- 
sure t o balance that of the ext e rnal medium (e.g. iMo rganti et al. 
19881; lHardcastle et alj Il998bl : IWorrall & Birkinshawl koOOh . In 
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pressure cannot be due solely to relativistic-electron dominance, as 
this would result in detectable X-ray inverse-Compton emission, 
which is not seen, while the presence of sufficient thermal gas at 
the temperature of the environment was also ruled out by the pres- 
ence of deficits in X-ray surface brightness at the positions of radio 
lobes. This is also the case for th e sample of cluste r-centre FR-I 
radio sources recently studied by iDunn et alj d2005l) . Relativistic 
protons intrinsic to the jets would require extremely large ratios of 
proton to electron energy in order to provide the missing pressure 
- such high ratios are implausible as they would lead to overpres- 
sured jets on kpc scales. The most plausible origin of the required 
additional pressure is either material that has been entrained and 
heated, or magnetic dominance of the lobes; however, these sce- 
narios are both difficult to test. 

Despite these uncertainties about radio-source energetics, our 
understanding of the impact of radio galaxies on their environ- 
ments has improved dramatically over the past few years. Direct 
evidence for heating by FR-I radi o galaxies has been found in a 
number of individual syst ems (e.g. lKraft et a l. 2003; Crost onet al.1 
|2007| ;I Forman et alj|2007l) . and statistical studies have shown that 
radio-source heating may have a significant impact in the oute r 
regions of galaxy groups dCroston et alj|200"5al ; Ijetha et alj|2007l) . 
Progress in hydrodynamical modelling has shown that a range of 
heating mechanisms may be capable o f heating the central regions 
of clusters to prevent cooling flows (e.g.lBasson & Alexandej2003l ; 
iBriiggen & Kaisej200ll ; lRevnolds et al j2002h and semi-analytical 
galaxy formation models that take into account radio-source heat- 
ing appear to be able to reproduce the properties of the nearby 
galaxy population more accurately than has previously been pos- 
sible dCroton et alj|2006l : iBower et al.| [2006). The growing accep- 
tance that AGN, and particularly low-power (FR-I) radio galaxies, 
play an essential role in galaxy and structure evolution means that it 
is crucial to understand which types of radio structures form under 
what environmental conditions, and how the environmental impact 
of these different structures varies. 

In this paper, we present new XMM-Newton observations of a 
sample of low-power (FR-I) radio galaxies, which range in mor- 
phology from tailed or plumed sources to those with lobes or 
bridges. We then combine this sample with previous observations 
and archive data to investigate the environmental properties of the 
nearby FR-I radio-galaxy population, to constrain their particle 
content and dynamics, and to study their environmental impact. 

Throughout the paper we use a cosmology in which Ho = 70 
km s _1 Mpc~\ fi m = 0.3 and Q.a = 0.7. Spectral indices a 
are defined in the sense S„ oc u~ a . Errors quoted on parameters 
derived from model fitting correspond to the la (68 per cent) con- 



fidence range for one parameter of interest, except where otherwise 
stated. 



2 DATA ANALYSIS 

We carried out XMM-Newton observations of a sample of five low- 
power (FR-I) radio galaxies, 3C296, 3C,76.1, NGC 1044, 3C31 
and NGC 315 t o follow up our earlier studies of 3C 66B, 3C449 
and NGC 6251 JCroston et alj2003l ; lEvans et al.l2005l) . In this sec- 
tion we describe the data analysis for the five new observations, 
plus one archive observation of NGC 4261 [the A GN emission de- 
tected i n this observation h as been described by Sambru na et al.l 
(2003); Gliozzi et al. ( 2003) and the extended emission was dis- 
cussed briefly in ICroston et al.1 d2005al) l. In Section [5] we com- 
bine these resu lts with our previou s a nalysis of 3C66B , 3C449, 
and NGC 625 1 ^Croston et al.l2003l and lEvans et al]2005t) . Tabled 
lists the observational details for the nine XMM-Newton observa- 
tions discussed in the present paper. For three of the radio galaxies, 
3C296, NGC 315 and NGC 6251 we also made use of previously 
published Chandra radial surface brightness profiles to improve our 
constraints on the inner gas distributions. Details of these profiles 
are given in Section[3] 

2.1 Radio data 

All of the radio galaxies in our sample apart from NGC 1044 are 
well-studied in the radio, with maps in the literature at multiple 
frequencies. Table [2] lists the radio maps used in the figures and 
analysis that follows. 

As existing low-frequency VLA observations of NGC 1044 
did not fully sample its large-scale structure, we observed it with 
the VLA at 20 cm to supplement existing archival data. In combi- 
nation with these archival data, the new observations were used to 
produce a good-quality map of of the large-scale radio structure for 
comparison with X-ray images. The VLA observations were car- 
ried out on 2004 March 19 in C configuration and 2004 July 17 in 
D configuration. The new data, together with an short archival ob- 
servation in the C configuration (AH276) were externally calibrated 
in the standard manner using AIPS. Each dataset was then imaged, 
cleaned and self-calibrated individually. The two C-configuration 
datasets were combined and imaged after first self-calibrating the 
archive dataset using a model from the new observations. The D- 
configuration data were then self-calibrated using the combined C- 
configuration model. Finally, the C- and D-configuration data were 
combined with appropriate weights and used to produce the images 
shown below. 



2.2 XMM-Newton data 

The data for the five new observations and one archival dataset not 
previously analysed were reprocessed using the XMM-Newton S AS 
version 6.0.0, and the latest calibration files from the XMM-Newton 
website. The pn data were filtered to include only single and double 
events (PATTERN < 4), and FLAG==0, and the MOS data were 
tilted according to the standard flag and pattern masks (PATTERN 
< 12 and #XMMEA_EM, excluding bad columns and rows). 

Several of the datasets were badly affected by background 
flares, and so GTI filtering was applied in order to ensure that ac- 
curate measurements could be made in low surface-brightness re- 
gions. Filtering was carried out by applying a ±20 per cent clip to a 
lightcurve in the energy range where the effective area for X-rays is 
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negligible (10-12 keV for MOS, 12-14 keV for pn). In one case, 
NGC315, initial analysis revealed significant contamination from 
soft proton flares, and so for this source a second stage of filtering 
was applied using the full energy range of 0.3 - 10 keV and an off- 
source region, so as to remove any additional soft flares, again by 
applying a ±20 per cent cut. 

To enable the use of double background subtraction tech- 
niques for spectral and spatial analysis (see Section [2~3l >, the events 
lists were vignetting corrected using the SAS task evigweight. 
Filter-wheel closed datasets (Pointecouteau et al., private commu- 
nication) for particle background subtraction were processed and 
filtered in the same way as the source datasets: the events lists for 
each camera were first filtered using the same FLAG and PAT- 
TERN filters as the data files, and weighted using evigweight. 
Weighting the particle background events lists allows this back- 
ground component to be subtracted correctly from the particle com- 
ponent in the source files, which must be weighted up in order 
that the vignetting correction is applied correctly to the source 
and background X-ray emission in the target observation. Appro- 
priate background files were then generated for each source using 
attcalc to recast the events list to the correct physical coordinates. 
The resulting background file appropriate to each source events file 
was used as part of the double subtraction process in the analysis 
described below. A scaling factor was determined for each back- 
ground dataset to account for differences in the normalisation of 
the particle and instrumental background between the source and 
filter-wheel closed datasets. The scaling factors were calculated by 
comparing the 10-12 keV (MOS) or 12-14 keV (pn) count rates for 
the source and background datasets. The background data products 
were scaled by this factor before carrying out background subtrac- 
tion of spectra or profiles. 

In the case of the pn camera, it is also necessary to take into 
account the contribution of out-of-time events: 6.3 percent of de- 
tected pn counts (in full-frame mode) are recorded with an incor- 
rect detector co-ordinate in the y direction. These events produce 
an additional position-dependent background component that can 
be important in the analysis of low surface brightness emission. We 
therefore generated an out-of-time events list for the pn camera us- 
ing epchain, which was scaled appropriately and subtracted as an 
additional background component in the analysis that follows. 



2.3 X-ray analysis 

2.3.1 Imaging and spectral analysis 

In order to make pictures of the extended emission and to make 
overlays on the radio emission we extracted images in the 0.5 - 5 
keV energy range from the three EPIC cameras for each dataset 
using evselect (these images were not used in our quantitative anal- 
ysis). We then combined the images to produce a single image for 
each source, weighting the two MOS images by a factor chosen 
to increase the MOS sensitivity to be equivalent to that of the pn 
data (this is necessary in order to remove chip-gap artefacts from 
the image). The combined image was then exposure-corrected us- 
ing an exposure map made by summing appropriate maps for each 
camera obtained using eexpm ap. This method is similar to that 
used by Bohring er et al.l J2007I) . The images were not corrected for 
vignetting as this excessively weights up the unvignetted particle 
background at large radii. Point sources were then identified and 
removed from the combined image and Gaussian smoothing was 
applied with a range of kernel sizes so as to emphasize structure 



on different size scales. X-ray images and radio overlays for each 
radio-galaxy environment are shown in Figs.[T]to|6] 

We carried out spectral analysis for several regions of each 
environment, as described in more detail in the following sections. 
Spectra were extracted using the SAS evselect task, with vignetting 
correction included by applying the event weights. We used a dou- 
ble subtraction technique that makes uses of filter-wheel closed ob- 
servations, which provide high signal-to-noise information about 
the particle and instrumental background as a function of spa- 
tial positiorQ. To account for variations in the level of particle 
background between our observations and the filter-wheel closed 
datasets, the background events were scaled by the ratio of 10-12 
keV (MOS) or 12-14 keV (pn) counts over the full field of view for 
the source and background datasets. Spectra were then extracted 
for source and local background extraction regions (typically a sur- 
rounding annulus) from both target and background events lists, 
with the background events scaled as explained above. The spec- 
tra extracted from the background events list were used as back- 
grounds for the corresponding target observation spectra. Appro- 
priate response files were generated using rmfgen and arfgen. The 
target spectra for both the source and local background regions 
were grouped to a minimum of 20 counts per bin after background 
subtraction. The local background spectrum from the target ob- 
servation was then used to determine the Galactic and cosmic X- 
ray background spectrum for our observation. By using the scaled 
filter-wheel closed spectrum for the same region as a background, 
the particle background is completely removed. We therefore fit- 
ted an X-ray background model to the local background spectrum 
consisting of two mekal models to account for emission from the 
Galactic bubble and a power-law model absorbed by the Galactic 
Nh in the direction of the target (see Table [6} to account for the 
cosmic X-ray background. The mekal temperatures wer e allowed to 
vary, but the power-law index was fixed at T — 1.41 jLumb et"al] 
2002). The normalisations of all three components were allowed 
to vary. In all cases this model provided a good fit to the back- 
ground spectrum of our target observations. For each source spec- 
trum, we used the corresponding filter-wheel closed spectrum as 
a background spectrum, to account for particle background, and a 
fixed X-ray background model consisting of the best-fitting model 
for the particular observation, with the normalisations of each com- 
ponent fixed at the best-fitting values scaled to the appropriate area 
for the source extraction region. This method gave results consis- 
tent with other methods (double background subtraction using the 
background template files of Read & Ponman 2003 and local back- 
ground subtraction), but in many cases the more accurate subtrac- 
tion of both particle and X-ray background components afforded by 
this method enabled us to achieve significantly better constrained 
temperatures and abundances. Spectral fitting results are summa- 
rized in Tableland discussed in detail in Section[3] 

2.3.2 Spatial analysis 

Analysis of the gas distributions in each object was carried out us- 
ing radial surface brightness profiles, extracted in concentric an- 
nuli with point sources and chip gaps masked out. In all cases 
where resolved jets have been detected with Chandra, they are 
within the core radius of the PSF so that they can be consid- 
ered to be part of the central point source, which is modelled as 

See the XMM-Newton background analysis website at 
http://xmm.vilspa.esa.es/external/xmm_sw_cal/background/index.shtml 



4 J. H. Croston et al. 

Table 1. Details of the new, archival and previously published XMM-Newton observations discussed in this paper. 



Object z Ang. scale Njj ObsID Date Exposure 

(kpc/arcsec) (10 20 cm~ 2 ) MOS 1 ,MOS 1 ,pn (s) 



NGC315 


0.016 


0.326 


5.88 


0305290201 


2005-07-02 


15236,14984,12580 


3C31 


0.017 


0.346 


5.39 


0305290101 


2005-08-03 


14734,14745,10819 


3C66B 


0.0215 


0.435 


8.36 


0002970201 


2002-02-05 


17915,17862,14692 1 


NGC 1044 


0.021 


0.425 


8.76 


0201860301 


2004-07-18 


13674,13839,13565 


3C76.1 


0.032 


0.639 


10.8 


0201860201 


2004-08-02 


17762,18165,10961 


NGC 4261 


0.0073 


0.150 


1.55 


0056340101 


2001-12-16 


26150,26191,16257 


3C296 


0.024 


0.484 


1.86 


0201860101 


2004-08-05 


19584,20315,14729 


NGC 6251 


0.0244 


0.492 


5.58 


0056340201 


2002-03-26 


23725,24715,8074 1 


3C449 


0.0171 


0.348 


11.8 


0002970101 


2001-12-09 


20865,20837,1670s 1 



1 The exposure times quoted for the three previously published datasets are l ivetimes after GT I-filtering using hard-band lightcurves, similar to the 
method used for the new observations, as described in ICroston et all i2003h and lEvans et"aD <2005h . 



Table 2. Details of the new and previously published radio maps used in the analysis presented in this paper. The radio map for NGC 426 1 was made from 
archival VLA data. 



Object 


Frequency (GHz) 


Observatory 


Date observed 


Reference/Proposal ID 


NGC 3 15 


1.365 


VLA 


2001 March 10 


Laing et al. (2006) 


3C31 


1.636 


VLA 


1995 April 28 


Laing et al. f2008) 


3C66B 


1.425 


VLA 


1994 November 10 


Hardcastle et al. (1996) 


NGC 1044 


1.425 


VLA 


2004 March 19 


AC712, AH276 


3C76.1 


1.477 


VLA 


1987 November 25 


Leahvetal. (1998) 


NGC 4261 


1.550 


VLA 


1984 April 23 


AP77 


3C296 


1.452 


VLA 


1987 November 25 


Leahy et al. (1998) 


NGC 6251 


0.337 


WSRT 


1988 May 12 


Leahvetal. f!998) 


3C449 


0.609 


WSRT 




Leahv et al. (1998) 



part of the surface brightness profile analysis. Vignetting correc- 
tion and double background subtraction using filter-wheel closed 
datasets was applied in the same manner as for the spectral analy- 
sis. The three surface brightness profiles were fitted separately with 
models including point-sour ce components, a single f3 model (e.g. 
iBirkinshaw & Worrall|[l993l) . and in some cases a model consist- 
ing of a line-of-sight projection of the sum of two (3 models in gas 
density, hereafter referred to as a projb model. The projb model is 
defined by: 



n(r) = n 



1 + 



-3/W2 



+ N 



r 

1 + — 

ri 



-3/3/2 



(1) 



where N is the relative normalisation of the two f3 model compo- 
nents. The corresponding surface brightness profile is determined 
via numerical integration of the following expression: 



S X (R) oc / n 2 (l,R)dl 



(2) 



where r — (I 2 + R 2 ) 2 is the radius, R is measured in the plane of 
the sky, and I is along the line of sight. We explored the three or six 
dimensional parameter space for the /3 and projb models respec- 
tively, using a Markov-Chain Monte Carlo (MCMC) methofl We 



2 Our implementation of the MCMC method uses the Metropolis- 
Hastin gs algorithm in a manner very similar to the METRO code of 



iHobson & Ba ldwin 1 2004), but is implemented to run on a cluster of multi- 
processor computers using the Message Passing Interface (MPI). 



used the joint \ value for the three profiles as the likelihood esti- 
mator and in the sections that follow we quote both the parameters 
of the best-fitting model and the Bayesian estimates for the value of 
each parameter. Plausible ranges for each parameter were estimated 
by carrying out extreme fits and these were used as priors for the 
MCMC method. The parameter ranges for each sou rce are given in 
Table [5] Uncertainties (formally credible intervals: Gregorv] [2005t) 
for each parameter are determined from the 1-dimensional projec- 
tion of the minimal n-dimensional volume that encloses 68 per cent 
of the posterior probability distribution as returned by the MCMC 
algorithm. The uncertainties determined in this way correspond to 
lcr errors for 2 or 5 interesting parameters for the j3 and projb 
models, respectively. Quantities derived from the surface bright- 
ness model fits, included X-ray luminosity and pressure at a given 
radius, are obtained by determining that quantity for each model fit 
and then obtaining the Bayesian estimate for the quantity in ques- 
tion. Uncertainties on the derived quantities come from the minimal 
one-dimensional interval that encloses 68 per cent of the posterior 
probability space for the given quantity, and so correspond to lcr 
errors for a single parameter of interest in a traditional fitting pro- 
cedure. As seen in Section [3] some of the parameters are strongly 
correlated, leading to large uncertainities on the individual model 
parameters; however, tight constraints can be obtained on derived 
quantities such as pressure and luminosity. Our method does not 
treat the model normalisation as a free parameter, instead finding 
the best normalisation for each model; however, the uncertainty in 
model normalisation does not significantly affect the constraints on 
model parameters or derived quantities. 
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Before fitting, each model was convolved with the 
XMM-Newton PSF based on the on-axis parametrization de- 
scribed in the XMM-Newton CCF files XRT1_XPSF_0006.CCF, 
XRT2_XPSF_0007.CCF and XRT3_XPSF_0007.CCF. For four 
sources with C/zanafra-detected galaxy-scale environments with 
core radii too small to be resolved by XMM-Newton (NGC315, 
3C296, and NGC6251), we included previously published 
Chandra profiles in the combined fits so as to model the inner 
profiles as well as possible. Final pn surface brightness profiles for 
each radio-galaxy environment with a model using the best-fitting 
model parameters, as discussed in the individual sections below, 
are shown in Fig. [7] 



3 RESULTS 
3.1 NGC315 

Images of the environment of NGC 315 are shown in Fig.Q] Emis- 
sion is detected from NGC315, the group galaxies NGC311 and 
NGC 3 18, a candidate more distant cluster at a = 00 ,l 57 m 31 s , 8 = 
+30° 13' , which does not appear to have been identified previously, 
and the intragroup medium. The X-ray emission appears somewhat 
elongated to the north-east in a direction perpendicular to the radio 
jets . We measure a total of ~ 4000 0.3 - 7.0 keV pn counts from 
the galaxy and group environment. Although NGC 315 has previ- 
ously been studied in detail with ROSAT and Chandra, this is the 
first detection of an environment on hundred-kpc scales. Unfortu- 
nately this dataset is also highly flare-contaminated, with only 12 - 
15ks usable of a total observed 50ks (see Table [6). 

The results of surface brightness profile fits for NGC 315 are 
in Table [4] The inner enviroment of NGC 315 is well constrained 
from Chandra observations, and so in this case we carried out a 
joint fit to four profiles : one for each XMM-Newton camera, plus 
the Chandra profile of IWorrall et al.l ( 120070 . As shown in Fig. [7] 
the profile of NGC 315 does not appear to flatten into a group-scale 
component on hundred kpc scales, unlike the majority of other FR- 
I environments known to date, and so the projb model was not 
required to fit its surface brightness profile. As shown in Table [4] 
the preferred model parameters are in agreement at the joint la 
level wi th the parame t er val ues measured from the Chandra data 
alone bv lWorrall et al.l J2007h (0 = 0.52±0.01 andr c = 1.7±0.2 
arcsec). 

We extracted a spectrum in the inner 60 arcsec to study the 
core emission from NGC 315 and found that the Chandra best- 
fitting model consisting of an absorbed power law with Nh = 
7.6 x 10 21 cm -2 and T = 1.57 plus a mekal model with kT = 
0.62 ± 0.02 keV and abundance fixed at 0.3 solar was a good fit 
(X 2 = 155 for 125 d.o.f.). The results of fitting to the group emis- 
sion in a region between 120 and 300 arcsec are given in Table[3] 

In order to obtain profiles of gas density and pressure we 
assumed a linear ramp in gas temperature obtained by fitting a 
straigh t line to the Chandra temperature profile of IWorrall et al.l 
(2007) and the XMM-Newton temperature data point. The tempera- 
ture was held constant at the XMM-Newton-measuved value beyond 
217 arcsec. Errors on profiles take into account both the uncertainty 
in model parameters and the the uncertainty in temperature. 

We were also able to extract spectra for the group member 
galaxy NGC 311. The spectra were well fitted by a mekal model 



between the two models; however, the temperature of the thermal 
fit is higher than would be expected for a galaxy atmosphere, and 
so it is likely that some of the emission is associated with nuclear 
activity or X-ray binaries. 



3.2 3C31 

Fig. [2] shows images of the environment of 3C31 from our new 
XMM-Newton observation. Emission is detected from several group 
galaxies as well as from the intragroup medium. Unfortunately, 
the observations of 3C 3 1 were severely fiare-contamined making 
quantitative analysis extremely difficult. After preliminary anal- 
ysis we concluded that it was not possible to estimate reliably 
the level of background for this observation and so we did not 
carry out either spectral or spatial analysis. 3C31 has been stud- 
ied in detail with ROSAT ( Komossa & Bohr ingeijl 19911 and Chan- 
dra dHardcastle et al.ll2002l) . For comparisons with the other radio 
galaxies in this sample, we therefore made use of the results of 
those two studies. A forthcoming reobservation with XMM-Newton 
will allow a much more detailed study of 3C 3 1 's environment and 
environmental interactions. 



with kT = 1.4" 1 



keV with abundance fixed at 0.3 times solar 



1.2 for 3 d.o.f.) or by a power law with T — 2.6 ± 0.5 



(X = 1-6 for 3 d.o.f.). The data are insufficient for distinguishing 



3.3 NGC 1044 

The environment of NGC 1044 is shown in Fig. [3] This dataset was 
relatively free of flare contamination, and a total of ~ 7840 0.3 - 
7.0 keV pn counts are detected from the environment. Fig. [3] re- 
veals a cavity at the position of the SE lobe. We tested the sta- 
tistical significance of the surface brightness deficit by compar- 
ing the total (MOS1 + MOS2 + pn) surface brightness in a cir- 
cle of radius 53 arcsec defined using the radio map, an annular 
background centred on the nucleus of NGC 1044 extending from 
the inner to the outer edge of the cavity region, but excluding 
the south-east quadrant, and a background of the same size ro- 
tated about the centre of the environment. The cavity has a surface 
brightness of 0.0247± 0.0017 counts arcsec -2 , compared to a sur- 
face brightness of 0.0315 ± 0.0005 counts arcsec -2 in the annular 
background region, and a surface brightness of 0.0398 ± 0.0022 
counts arcsec -2 in the same-sized background region. The errors 
on surface brightness measurements were obtained by combining 
in quadrature source and background errors determined from the 
Poissonian errors on the counts in each region. The deficit is signif- 
icant at a level > 3a. 

NGC 1044 has not previously been observed in the X-ray, 
and so its inner (galaxy-scale) environment is not well-constrained. 
Our new observations provide the first detection of an X-ray envi- 
ronment for this radio galaxy. The surface brightness profile for 
NGC 1044 (Fig. [7} shows a characteristic two-component shape 
indicating the presence of a central AGN or galaxy-scale compo- 
nent as well as a large-scale group gas component. A point-source 
plus single model is a poor fit to the profiles, and so we used a 
point-source plus projb model, with model fitting results given in 
Table [4] Results of a spectral fit to the group-scale emission in an 
annulus between 100 and 400 arcsec are given in Table [3] 

Although the surface brightness profile does not show evi- 
dence for a dominant AGN component, we extracted spectra from a 
source-centred circle of one arcminute radius to investigate the con- 
tributions of AGN and galaxy-scale gas-related emission. Consis- 
tent with the surface brightness profile results we found that a sin- 
gle power-law model did not provide an acceptable fit to the data. 
We obtained a good fit (x 2 = 40 for 37 d.o.f.) for a single mekal 
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Figure 1. The X-ray environment of NGC315. Images made from the combined MOS1, MOS2 and pn events lists in the 0.5 - 5.0 keV energy range, with 
exposure correction to correct for chip gaps but not vignett ing, as described in the text. Left: smoothed with a Gaussian of FWHM 33 arcsec to show the 
group-scale emission. Right: with 1.4-GHz radio contours (Laing et al. 2006) overlaid. 
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Figure 2. The X-ray environment of 3C 31. Images made from the combined MOS1, MOS2 and pn events lists in the 0.5 - 5.0 keV energy range, with exposure 
correction to correct for chip gaps but not vignetting, as described in the text. Left: smoothed with a Gaussian of FWHM 1 3 arcsec to highlight the emission 
associated with several of the NGC 383 group galaxies. The central galaxy is NGC 383, with a slight elongation to the south west due to emission from its 
companion galaxy NGC 382; right: s moothed with a Gau ssian of FWHM 26 arcsec to highlight the group-scale emission, and with galaxy emission removed, 
overlaid with 1.6-GHz radio contours ILaing et al . 2008). 



model with kT = 0.82 ± 0.05 keV and Z = O.15I^ 5 Z . A 
somewhat better fit was obtained for a two-component mekal plus 
power-law model (x 2 = 32.3 for 36 d.o.f.) with kT = 0.77t° ° l 
keV, abundance fixed at 0.5 times solar and T — 1.91^ ^. The 
single mekal model gives an unabsorbed 0.1 - 10 keV luminos- 
ity of (1.3 ± 0.2) x 10 41 ergs s -1 , which is high for a galaxy 
atmosphere, whereas we obtain a more reasonable luminosity of 
7.5 x 10 40 erg s" 1 for the combined fit. The combined fits gives 
a 1-keV flux density of 6 ± 1 njy for the power-law component, 
which is cons istent with the nuc lear X-ray emission predicted by 
the relation o f lEvans et alj d2006l) for a source of NGC 1044's core 
radio flux density (< 0.08 Jy at 1.5 GHz). 

We used the surface brightness profile models and the mea- 



sured temperature to determine profiles of gas density and pressure 
by using the best-fitting projb model parameters to obtain an emis- 
sion measure profile. As we were unable to obtain a temperature 
profile for this group, we assumed the gas is isothermal with the 
spectral parameters of the best fit to the group spectrum (kT = 1.3 
keV and Z — 0.12 solar) in order to obtain gas density and pres- 
sure profiles. Errors on the density and pressure profiles were de- 
termined by combining in quadrature the la confidence range from 
the dispersion in profile fits with the uncertainty in temperature as 
for the other sources. 

Of the NGC 1044 group galaxies, X-ray emission is clearly 
detected from NGC 1046 (z = 0.0199) and CGCG 414-036 
(z = 0.0206), as shown in Fig. [3] We were able to extract spectra 
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Figure 3. The X-ray environment of NGC 1044. Images made from the combined MOS1, MOS2 and pn events lists in the 0.5 - 5.0 keV energy range, with 
exposure correction to correct for chip gaps but not vignetting, as described in the text. Top left: smoothed with a Gaussian of FWHM 13 arcsec to highlight 
the emission associated with several of the NGC 1044 group galaxies; top right: the same image smoothed with a Gaussian of FWHM 20 arcsec and with 
galaxy emission removed to show the group-scale emission and cavity to the S; an inset shows a close-up of the cavity region from an image smoothed with 
a Gaussian of FWHM 26 arcsec; bottom left: image smoothed with a Gaussian of FWHM 26 arcsec to highlight the group emission, with 1.4-GHz radio 
contours overlaid illustrating that the southern radio plume fills the cavity shown in the middle panel; bottom right: regions used to investigate significance of 
cavity detection. 



for NGC 1046, and fitted mekal, single power-law and mekal plus 
power-law models. All three models gave good fits to the spectra; 
however, the best-fitting single mekal model (x 2 = 4.2 for 6 d.o.f.) 
required low abundance (Z < Q.QAZq) and the best-fitting power- 
law model (x 2 = 1-8 for 7 d.o.f.) required a steep photon index 
(T = 3.2^0 *). We obtained the lowest reduced x 2 an d most phys- 
ically plausible parameters for the two-component fit (x 2 = 0.8 
for 5 d.o.f.) with kT = 0.2±£? keV and T = 2.3±J;? and abun- 
dance fixed at 0.5 times solar. We conclude that there is likely to be 
some non-thermal emission associated with this galaxy, either asso- 
ciated with X-ray binaries or nuclear activity, in addition to thermal 
emission from a galaxy atmosphere. There were insufficient counts 
associated with CGCG414— 036 to extract a useful spectrum. 



3.4 3C76.1 

This observation is the first pointed X-ray observation of 3C76.1. 
As shown in Fig. [4] the environment of 3C76.1 is considerably 
poorer than those of the other radio galaxies in this sample. The 
X-ray emission appears elongated in the north-south direction, per- 
pendicular to the radio lobes. We detect a total of ~ 385 0.3 - 7.0 
keV pn counts from the group environment. 

To confirm the presence of extended emission, we fitted a 
point-source-only model to this surface brightness profile (which 
is shown in Fig [7}. This did not give an acceptable x 2 , whereas 
we were able to obtain good fits with a point-source plus f3 model, 
as listed in Table [4] As the profiles are inconsistent with a point- 
source only model, and the /3-model fit is good, we conclude that 
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3C76.1's environment has been firmly detected for the first time. 
As discussed below, the point source component appears to orig- 
inate partly from a central AGN and partly from an unresolved 
galaxy-scale atmosphere. 

We list the results of a spectral fit to the group emission in the 
region between 30 and 200 arcsec in Table[3] We used a fixed abun- 
dance as the abundance value tended to unrealistically high values 
if left free. We used the fitted surface brightness profile model and 
the measured temperature to determine profiles of gas density and 
pressure. As for NGC 1044, we assumed the gas to be isothermal 
with the spectral parameters of the best fit to the group spectrum 
(kT — 0.91 keV and Z = 0.3 solar) in order to obtain gas density 
and pressure profiles. Errors on the profiles were determined as for 
the other sources described above. 

To constrain the origin of the unresolved central emission, 
we extracted a spectrum in the central 30 arcsec using local back- 
ground to subtract off the contribution from the group-scale emis- 
sion. We found that neither a power law nor a mekal model on 
their own gave an acceptable fit, and so we fitted a power-law plus 
mekal model, which gave a \ 2 value of 15.5 for 12 d.o.f. with 
kT = 0.43±S:1? keV ( for abundance fixed at 0.3 times solar) and 
r = 0.9lg'g. Such a flat power-law index suggests that an ab- 
sorbed power-l aw component, as commonly seen in narrow-line 
radio galaxies dHardcasfle. Evans & Crosto n 2006) a nd in some 
low-power sources such as NGC 3801 dCroston et al.ll2007l) . may 
be present. 

3.5 NGC 4261 

The environment of NGC 4261 is shown in Fig. [5] The inner en- 
vironment of_jhis_£adiog_alaxy a s observed by Chand ra was dis- 
cussed bv lZezas et alj d2005t) and lCroston et al.l d2005al) . where we 
also presented a n image of the XM M-Newton data but not a detailed 
analysis, and bv lJefha et ail d2007t) . We detect a total of ~ 33, 000 
0.3 - 7.0 keV pn counts from the source, including a contribution 
from the central AGN. The dataset is not significantly contaminated 
by flares. 

Fig. [5] reveals a cavity at the position of the W lobe. We 
tested the statistical significance of the surface brightness deficit 
by comparing the total (MOS 1 + MOS2 + pn) surface brightness 
in a circle of radius 67 arcsec defined using the radio map with 
an annular background region centred on the nucleus of NGC 4261 
extending from the inner to the outer edge of the cavity region, 
but excluding the south-east quadrant, and with a background of 
the same size as the cavity region but rotated about the centre 
of the group to the south. The cavity has a surface brightness of 
0.0942±0.0027 counts arcsec -2 , compared to a surface brightness 
of 0.1081 ± 0.0011 counts arcsec -2 in the annular background re- 
gion, and a surface brightness of 0.1434 ± 0.003 counts arcsec -2 
in the same-sized background. Hence the deficit is significant at a 
level between 3 and 8<j. The errors on the surface brightness mea- 
surements were obtained as for NGC 1044. There is no statistically 
significant large-scale cavity on the eastern side of the source; how- 
ever, there is evidence for a smaller region of lower surface bright- 
ness deficit close to the nucleus and coincident with the inner parts 
of the eastern radio lobe. 

As shown in Fig. [7] the surface brightness profile shows a 
two-component form with a large-scal e component in a ddition to a 
galaxy-scale inner environment (e.g. IZezas et al]|2005h . We there- 
fore fitted the point-source plus projb model as for NGC 1044. The 
results are listed in Table|4] Table[3]gives the results of spectral fit- 
ting to the group-scale emission between 60 and 600 arcsec. We 



were concerned that because the physical scales we are sampling 
in NGC 4261 are considerably smaller than in the other sources in 
our sample (due to its lower redshift) contamination from X-ray 
binaries in the galaxy could be a problem (D25 = 4.1 arcmin 
for NGC 4261, whereas it is < 2 arcmin for most of the other 
galaxies - here D25 is the optical radius of the galaxy at a B- 
band surface brightness of 25 mag arcsec 2 ), and so we included 
a thermal bremsstrahlung component with kT — 5.0 keV to ac- 
count for an X-ray binary population; however, for the global spec- 
trum this component was not significant. We investigated whether 
a temperature gradient is present by extracting spectra in several 
annuli. For the outer two bins the normalization of the included 
thermal bremsstrahlung component was negligible, consistent with 
the expectation of no contribution from binaries on those scales. 
In the inner two regions (which are both within the D25 radius of 
the galaxy) the bremsstrahlung component contained a significant 
fraction of the total flux, corresponding to a total X-ray binary 0.3 
- 8 keV luminosity of ~ 3 x 10 40 erg s -1 . This is consistent with 
the expectation for NGC 4261 based on the mean Lx /Lb relation 
of Kim & Fabbiano (1994). We found no evidence of a signifi- 
cant temperature gradient when the X-ray binary contribution is 
accounted for. 

We also extracted spectra from a source-centred circle of one 
arcminute radius in order to investigate the AGN and galaxy-scale 
X-ray emission. We found that neither a single mekal nor a sin- 
gle power-law model provided an acceptable fit. We then fitted a 
mekal plus power-law model, which gave a reasonable fit statistic 
(X 2 = 503 for 408 d.o.f.) for kT = 0.65±g;^ keV and abun- 
dance fixed at 0.5 times solar; however, the best-fitting power-law 
index is very flat: V = 0.57 ± . 11. W e therefore tried fitting the 
best-fitting model of IZezas et ai] d2005l) (their Model 5). If we fix 
the column densities, power law indices and Raymond-Smith tem- 
perature to their best-fitting values, we cannot obtain an acceptable 
fit; however, allowing the temperature to vary gives a good fit with 
X 2 = 585 for 406 d.o.f. for a temperature of kT = 0.69 ± 0.01 
keV. The higher temperature is likely to be due to the larger ex- 
traction region used here, which is likely to include hotter gas from 
outer regions of the galaxy atmosphere and/or a contribution from 
X-ray binaries. 

Two other galaxies in the NGC 4261 group were detected 
with sufficient counts to carry out spectral fitting, NGC 4254 (2 = 
0.016) and PCG 039660 {z = 0.017). We fitted single mekal and 
power-law fits to each galaxy (the statistics were too poor to allow 
fitting of more complex models). For NGC 4254 we found good 
fits with both models, with a slightly better fit statistic (x 2 = 5.9 
for 4 d.o.f.) for the mekal model, but a temperature of > 1.3 keV, 
which is high for a galaxy atmosphere. The best-fitting power-law 
model (x 2 = 6.6 for 4 d.o.f.) had V = 1.4+{J;Jj. It seems likely that 
the emission from this galaxy is at least partly non-thermal in ori- 
gin. For PCG 039660, the mekal model was unacceptable, whereas 
the power-law model was a good fit (x 2 = 2.0 for 2 d.o.f.) with 
r = 1.61q' 6 , and so we conclude that its emission is predomi- 
nantly non-thermal. 

3.6 3C296 

As shown in the left hand panel of Fig. [6] 3C 296 possesses a group 
environment of moderate X-ray surface brightness. This is the first 
detection of a group-scale environment for this radio galaxy. We 
measure a total of ~ 6800 0.3 - 7.0 keV pn counts from the group 
environment. 

Fig.|6]reveals a cavity at the position of the NE lobe. We tested 
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Figure 4. The X-ray environment of 3C76.1. Images made from the combined MOS1, MOS2 and pn events lists in the 0.5 - 5.0 keV energy range, with 
exposure correction to correct for chip gaps but not vignetting, as described in the text. Left: image smoothed with a Gaussian of FWHM 5.2 arcsec; right: 
smoothed wih a Gaussian of FWHM 10.4 arcsec and with 1.4-GHz radio contours overlaid. 



the statistical significance of the surface brightness deficit by com- 
paring the total (MOS1 + MOS2 + pn) surface brightness in a circle 
of radius 74 arcsec defined using the radio map, in an annular back- 
ground region centred on the nucleus of 3C 296 extending from the 
inner to the outer edge of the cavity region, but excluding the north- 
east quadrant, and in a same-sized circular region rotated about the 
centre of the group to the south. The cavity has a surface brightness 
of 0.0736 ± 0.002 counts arcsec -2 , compared to a surface bright- 
ness of 0.0870 ± 0.001 counts arcsec - in the annular background 
region, and a surface brightness of 0.0927± 0.002 counts arcsec -2 
in the same-sized region. The errors on the surface brightness mea- 
surements were obtained as for NGC 1044. Hence the deficit is sig- 
nificant at a level of ~ 4.5a. 

3C 296 has pre viously been observed with Chandra 
dHardcastle e7ai1l2005l) . These observations revealed the presence 
of a galaxy-scale hot-gas component with kT = 0.7 keV ex- 
tending to distances of ~ 1 arcmin. This component dominated 
over the compact emission associated with the AGN beyond ~ 1 
arcsec. Our XMM-Newton surface brightness profile for 3C296 
(top lefthand panel of Fig. [7} is dominated by the galaxy-scale 
component to a distance of ~ 60 arcsec, beyond which the profile 
flattens slightly. 

As for NGC 315, we carried out a joint fit to the Chandra and 
XMM-Newton profiles for 3C 296, so as to provide the best possi- 
ble constraints on both the inner and outer gas properties. A single 
B model is an unacceptable fit to the profiles, due to the excess 
of emission at large radii. We therefore used the point-source plus 
projb model to model both the inner and outer components, find- 
ing good fits, as listed in Table|4] The inner parameters we measure 
are in good agreement with the best-fitting parameters of a single 
B model fitted to the Chandra profile alone, B = 0.56 ± 0.01 and 
r c = 0.7+o jo Mcsec dHardcastle et alj|20050 . Table[3]lists the re- 
sults of a spectral fit to the group emission in the region between 
50 and 400 arcsec. As the fit is comparatively poor, we extracted 
spectra in concentric annuli to investigate whether a temperature 
gradient is present. The temperature was found to decrease from 
~ 1.4 keV towards the centre to ~ 0.9 keV in the outer regions; 
however, the temperatures are consistent within the errors, so that 



no significant gradient can be claimed. We used the fitted surface 
brightness profile model to determine profiles of gas density and 
pressure, assuming an isothermal atmosphere. Density and pressure 
profiles, and their uncertainties were determined as for the other 
sources. While the Chandra observation of 3C 296 showed a signif- 
icant temperature gradient in the inner regions of the galaxy-scale 
emission, the inner temperature beyond 4 arcsec differs by at most 
10 per cent from the XMM-Newton global temperature, and so this 
difference is accounted for by the comparatively large uncertainty 
of the XMM-Newton global temperature, which we include in the 
pressure profile uncertainties. 

We also extracted a spectrum for the core region and for the 
nearby group galaxy NGC 5 531 (z = 0.026), shown to the north- 
west in Fig. [6] As found by lHardcastle et al.1 d2005h a power-law 
with only Galactic absorption was not a good fit to the core spectra. 
In good agreement with the Chandra results, we found a good fit 
to an absorbed power-law plus mekal model, with kT — 0.74l '5 
keV, N H = (2± 1) x 10 21 cm -2 , and T = 1.6±0.4. The gas tem- 
perature is somewhat higher than that measured by Chandra, which 
is not surprising as our region encompasses more of the galaxy at- 
mosphere, which is hotter on scales of tens of kpc. For NGC 553 1 , 
we found a good fit to a thermal model, with kT — 0.60+o'^ 2 keV 
and Z = O.2± U ;£Z0 {\ 2 = 17.2 for 14 d.o.f.), in contrast to the 
results with Chandra, which required the inclusion of the power- 
law component. When a power-law component was included, its 
normalization tended to negligible values, and the fit was not im- 
proved. We measure an unabsorbed 2- 10 keV flux density a factor 
of ~ 17 lower than that reported by Hardcastl e et alj (2005), and 
so variability of an AGN component may be the explanation for 
the differing spectral results. 

We can use the observed cavity associated with the eastern 
lobe to obtain an estimate for the jet power for 3C 296. If we use 
the sound crossing time to the edge of the lobe (2 x 10 8 yr) as an 
upper limit on the source age, and calculate the energy required to 
inflate the lobe as 4:P m idV, where P m id is the external pressure 
at the midpoint of the lobe (see Table[6}, and V is the volume es- 
timated from the extent of radio-lobe emission (this is preferable 
to estimating the cavity size directly from the X-ray data, due to 
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Figure 5. The X-ray environment of NGC 4261. Images made from the combined MOS1, MOS2 and pn events lists in the 0.5 - 5.0 keV energy range, with 
exposure correction to correct for chip gaps but not vignetting, as described in the text. Top left: smoothed with a Gaussian of FWHM 26 arcsec and with galaxy 
emission removed to show the group-scale emission and cavity to the W; top right: smoothed with a Gaussian of 33 arcsec and with 1.4-GHz radio contours 
overlaid illustrating that the W lobe coincides with the cavity in the group-scale gas; bottom: regions used to investigate significance of cavity detection. 



XMM-Newton's spatial resolution and the need for smoothing of 
the X-ray data to reveal the cavity structure), we find that the en- 
ergy required to innate the lobe is ~ 3.9 x 10 51 J (3.9 x 10 58 ergs), 
corresponding to a jet power of 1.2 x 10 36 W (1.2 x 10 43 erg s _1 ) 
for the two jets, assuming the western jet has the same power. This 
inferred jet power is roughly an order of magnitude higher than the 
X-ray luminosity of the 3C 296 group, so that significant heating 
is occuring in addition to the energy input required to balance ra- 
diative lo sses. The cavity properties ar e within the range shown in 
Fig. 8 of iMcNamara & Nulseiil < l2007t) , but lie to the high cavity 
power side of the cavity power/X-ray luminosity distribution. This 
jet power estimate is an order of magnitu de lower than the estimat e 
for 3C31 from the kinematic model of lLaing & Bridle! {2002a), 
which is somewhat surprising given that 3C 296 has a higher lu- 
minosity at 178 MHz (see Table |6). The sound crossing time is 
likely to be an overestimate of the lobe expansion time, as the lobes 
appear to be slightly supersonic in their outer regions at the present 
time, and were probably more highly overpressured at earlier times. 



Hence our calculation of the energy input from 3C 296 discussed 
above may be an underestimate of the true value. 



4 RESULTS FOR 3C 66B, 3C 449 AND NGC 6251 

For the three radio galaxies where we have previously published 
detailed analyses of the environmental properties as probed by 
XMM-Newton we ch ose to make use of the results published in 
ICroston et al. | te003h and lEvans et alj [2005), where very similar 
analysis procedures were used. In Table [3] we list the global gas 
properties obtain ed from single therm al model fits to the group 
gas as reported in lCroston et alj (2003) and Evan s et alj J2005h . To 
enable a fair comparison with the pressure profiles for the newly 
observed objects, we refitted the surface brightness profiles for 
the three previously observed sources with the point-source plus 
projb model where an inner component is detected (3C 66B and 
NGC 6251) or a point-source plus f3 model (3C449) and used the 
Markov-Chain Monte Carlo method, which improves the sampling 
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Figure 6. The X-ray environment of 3C296. Images made from the combined MOS1, MOS2 and pn events lists in the 0.5 - 5.0 keV energy range, with 
exposure correction to correct for chip gaps but not vignetting, as described in the text. Images are smoothed with a Gaussian of FWHM 19 arcsec. Top left: 
group emission; top right: cavity to the NE highlighted; bottom left: with 20-cm radio contours overlaid illustrating that the northern lobe fills the cavity seen 
in the top right panel; bottom right: regions used to investigate significance of cavity detection. 



of par ameter space . For N GC 6251, we included the Chandra pro- 
file of lEvans et all te005h to better constrain the small core-radius 
inner component. The results of these fits are listed in Table [4] We 
note that the quality of the fits for 3C 66B and 3C 449 is poor, as 
w as the case for the sin gle f3 model fits to the outer regions reported 
in lCroston et al.l {2003). We attribute this to the very high signal-to- 
noise data for these two sources, so that real small-scale deviations 
in the gas distributions may be affecting the quality of the fit. 



5 DISCUSSION 

Table [6] summarizes the group properties as measured from the 
XMM-Newton data. The three previously studied radi o-galaxy en- 
vironm ents are included, wit h properties taken from lCroston et all 
d2003h and lEvans etaflfcOOSl) . Properties of the radio galaxies ob- 
tained from the radio maps and from the literature are also included 
for comparison. It is clear from these results and from Table[3](and 
can also be seen from the X-ray images) that the low-power ra- 



dio galaxies in this "representative" sample inhabit a wide range of 
environments, with bolometric X-ray luminosities spanning nearly 
two orders of magnitude fr om ~ 10 41 erg s -1 to 10 43 e rg s _1 , con- 
sistent with earlier results jWorrall & Birkin shaw 2000). The result 
is important for theoretical models of radio-galaxy evolution and 
environmental impact, as it has sometimes been assumed by mod- 
ellers that environmental properties do not vary significantly. 

In the following sections, we discuss the relationship between 
environmental properties and radio structure, source dynamics and 
particle content. 

5.1 Relationship between properties of the radio galaxy and 
of the environment 

The importance of jet/environment interactions for group and clus- 
ter gas has been highlighted by the large number of X-ray cav- 
ity systems now kno wn, ranging from el l iptical galaxies to the 
richest clusters (e . g. iBohringer et aHll993l; iDunn & Fabianll2004l ; 
iBirzan et al]|2004 iMcNamara & Nulsenll2007t) . Of the nine sys- 
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Figure 7. EPIC pn surface brightness profile with best-fitting point source plus (3 model or projb model, as discussed in the text, for the newly observed 
galaxies. Top row (1-r): the bridged sources 3C 296, 3C 76.1, NGC4261; bottom row: the plumed sources NGC 315 and NGC 1044. (The somewhat poor fit in 
the outer regions for NGC 4261 is because the model is from the best joint fit, rather than the best fit to the pn profile alone). 



Table 3. Spectral fits to the large-scale group emission for each environment. Results are for joint fits to the three XMM-Newton cameras in the energy range 
0.3 - 7.0 keV. Luminosities are the Bayesian estimates for the bolometric luminosities to a radius of r'500 (see Table[6) with lcr errors, as determined from our 
surface brightness profile model fits. 
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Notes 


NGC 3 15 


0.016 


5.88 


1 0+0.5 


0.3 (fixed) 


117(84) 


41.89 + 0.01 




3C31 


0.017 


5.39 


1.5 ±0.1 


0.35 (fixed) 


25 (22) 


43.10 + 0.01 


Results from OP04 


3C66B 


0.0215 
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1 7 ,+0.03 


0.28 + 0.03 


819 (718) 


43.03+. J^ 3 


See C03 for details 


NGC 1044 


0.021 


8.76 


1 13+0.14 


fi4 +1 ' 3 


375 (330) 


42.13 + 0.08 




3C76.1 


0.032 


10.8 


oql +0.25 


0.3 (fixed) 


26.6 (25) 


41 30+002 




NGC 4261 


0.0073 


1.55 


1 45+0-23 
^^-O.Ul 


26 +0 15 


1186 (819) 


42 n7 +0 ' 21 




3C 296 


0.024 


1.86 


90+° 10 


14+ 006 


1071 (728) 


42 40+ - 13 




NGC 6251 


0.0244 


5.82 


1 6 +0 - 5 


4 +0 ' 7 


119(109) 


41 03+0.02 


See E05 for details 


3C449 


0.0171 


11.8 


0.98 + 0.02 


13+ 001 
u - lo -0.02 


881 (441) 


43.08 ± 0.04 


See C03 for details 



References: OP04 refers to Osmond & Ponmanl E(xA : C03 refers to Croston et ai] j2003l) : E05 refers to Eva ns et alj ^20051) . 



terns studied here with XMM-Newton, 6 have significant detections 
of cavities associated with the radio lob es [3C 296, NGC 1044, 
NGC 4261, 3C66B dCroston et al] l2003h . 3C449 JCroston et all 
120031) and NGC 625 l lE vans et alj|2005l) l. In several cases where 
a cavity is not identifiable in the X-ray data, the emission is nev- 
ertheless elongated perpendicular to the direction of the radio axis 
(3C 76. 1, NGC 315), which suggests a link between the structure of 
the X-ray emitting gas and the radio source. 



5.1.1 Correlations between X-ray and radio properties 

While the images presented in the previous sections clearly indicate 
an important relationship between the radio and X-ray properties 
of these systems, the extent to which the environmental properties 
control the large-scale radio morphology remains difficult to estab- 
lish observationally. We looked for correlations between the radio 
luminosity, radio-source projected physical size, axial ratio, inter- 
nal pressure, and the environmental properties as determined from 
our X-ray observations (X-ray luminosity, gas temperature, /3, r c , 
external pressure, ratio of external to internal pressure, and pres- 
sure gradient at the radius of the lobes). Table [7] lists the results 



FR-I environments with XMM-Newton 



13 



10 100 
Radius (kpc) 



1000 HD.1 



3C 66B 




Radius (kpc) 



10 100 
Radius (kpc) 



1000 0.1 



10 100 
Radius (kpc) 



NGC 1044 


a 
'a 

Q_ « 
o b 


NGC 6251 ; 










& a 










\ 

both \ 


o 








: w - 




'o 





10 100 
Radius (kpc) 



Radius (kpc) 



Figure 8. Pressure profiles for the XMM-observed radio-galaxy environments. The grey shaded area indicates the la uncertainty, which takes into account 
uncertainties in all the model parameters and in the temperature measurement. Hatched regions indicate an extrapolation of the pressure profile model beyond 
the region covered by the data. Horizontal lines indicate internal pressures of radio lobes. Top row (1-r): the bridged sources 3C296, 3C76.1, NGC 4261; 
middle row: the plumed sources 3C66B, NGC 315, and 3C449; bottom row: NGC 1044, and NGC 6251. 



of Spearman rank tests for each pair of parameters. The strongest 
correlations are plotted in Fig. 1101 

We found that radio luminosity is uncorrelated with any of the 
group gas properties; however, some correlations were found be- 
tween external pressure (measured at the midpoint of the source) 
and the projected radio-source size, axial ratio and external pres- 
sure. In particular, correlations that are significant at the ~ 90 per 
cent confidence level or higher were found between radio source 
size (in kpc) and both external pressure and pressure gradient (null 
hypothesis probabilities of 1.2 per cent and 7 per cent, respec- 
tively), and between axial ratio and both external pressure and pres- 
sure gradient (null hypothesis probabilities of 5.6 per cent and 9.1 
per cent, respectively). Although these results are not highly sig- 
nificant (and are somewhat dominated by the two largest sources, 
NGC 315 and NGC 6251: see Fig.HOt. they suggest a link between 



radio-source evolution and group gas properties. However, the cor- 
relation of radio-source size with external pressure and pressure 
gradient may simply result from a selection effect: we might expect 
the largest sources to probe further out into the group atmospheres, 
where the pressure is lower and the gradient likely to be steeper. To 
test this explanation, we determined the external pressures at a fixed 
physical radius of 300 kpc (comparable to the projected distances 
probed by the lobes of NGC315 and NGC 6251), which we con- 
sider to be an effective pressure normalization for the outer regions 
of each group. If the correlation between current source size and 
axial ratio with external pressure acting on the lobe is due solely to 
the fact that the larger sources can reach lower pressure regions of 
the groups, then we would expect no systematic difference in the 
pressure normalization at a fixed physical radius as a function of 
source size or axial ratio. Indeed, this is what we find, and so we 
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Table 4. Results of surface brightness profile fits. For each parameter we list the parameter value from the fit with the minimum \ 2 and the Bayesian parameter 
estimate together with credible intervals taking into account all fit parameters. Column 2 gives the model normalizations in electron density for the best fit in 
units of 10 3 m~ 3 - note that for the projb model this does not correspond exactly to the density at zero radius. Radii are given in units of arcsec. 



Object 



Best fit 







P 


r c 




^c, in 


N 


X 2 (dof) 




NGC315 


400 


0.54 


2.12 








186(126) 


0.54 


3C66B 


12 


1.2 


428.2 


1.2 


17.6 


0.083 


470 (230) 


0.91 


NGC 1044 


14 


1.2 


416 


0.84 


13.3 


0.0225 


38.2 (28) 


0.81 


3C76.1 


2.6 


0.94 


42.1 








11.3 (16) 




NGC 4261 


77 


0.30 


243.6 


0.64 


7.0 


0.00443 


156 (100) 


1 


3C296 


720 


0.31 


77.5 


0.60 


0.95 


0.0000446 


80.5 (92) 


i 


NGC 6251 


38 


1.1 


150 


1.2 


9.9 


0.044 


160 (139) 


l.( 


3C449 


3.7 


0.42 


57.1 








457 (273) 


0.42 



Bayesian estimate 

Pin 



N 



+0.02 
-0.01 
+0.29 
-0.47 
+0.39 
-0.51 

0.9* 

0.62* 

0.67* 
+0.2 

+0.06 
-0.05 



2.13 
367 
356 



+0.45 
0.39 
-103 
-155 
+640 
-255 
43* 

482* 
545* 

58 ii6 



1.1 

0.88 



+0.1 
0.4 
+0.32 

0.35 



64+ 35 

U ' OB -0.03 
0.96* 



16.4 



+3.6 
10.1 



14.2"* 



70 +3.0 

sfi +0.44 
■ 8b -0.31 

8.1* 



0.094 
0.0243 



0.3 
0.07 
+0.0254 
0.0185 



0.0046* 

00091+ - 0008 
u.uuuzi_ 0001 

0.053 ± 0.050 



Unconstrained within the parameter range explored. See Table|5] 



Table 5. Bayesian priors for the surface brightness profile fits. Priors in core radius and relative normalisation were determined individually for each source 
by finding extreme fits that gave poor \ 2 values. 



0min Pmax ^c,min ^c,max 0in,min Pin,max ^cin,mm ^cin.max ^min N ma x 



NGC 315 
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1 


30 
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0.5 
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50 
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0.3 


1.2 
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20.0 
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0.02 
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1.2 


2 


100 
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0.3 


1.2 


50.0 


2000 


0.3 


1.2 


0.1 


10.0 


0.00001 


0.001 


3C296 


0.3 


1.2 


50 


1000 


0.3 


1.2 


0.1 


10.0 


0.00001 


0.001 


NGC 6251 


0.3 


1.2 


15 


500 


0.3 


1.2 


15 


500 


0.0001 


0.1 


3C449 


0.3 


1.2 


10 


200 















conclude that the correlations with external pressure and its gradi- 
ent are most likely to be a selection effect. The lack of correlations 
between radio source structure and f3 and r c support this explana- 
tion. As radio source size and axial ratio are themselves strongly 
correlated (larger sources tend to have small axial ratios), it is plau- 
sible that the correlations with axial ratio have the same explana- 
tion; however, the size-axial ratio relationship itself is likely to be 
telling us something about FRI evolution, and so the different ex- 
ternal pressures and pressure gradients experienced by the smaller, 
and typically fatter sources (high pressures and flatter pressure gra- 
dients) compared to the larger, and typically thinner sources (low 
pressures and steeper pressure gradients), may be a determining 
factor in evolving a large FRI source. 

We also find significant correlations between both radio- 
source size and axial ratio and the pressure ratio, P ex t/Pint, with 
null hypothesis probabilities of 2.9 per cent and 10 per cent, re- 
spectively. We will discuss the relationship between pressure ratio 
(which we consider to be a measure of the size of the energetic 
contribution required from non-radiating particles - see below) and 
radio-source structure in more detail in Section [5"!2l 

The lack of correlations between X-ray luminos- 
ity/temperature and any of the radio properties of these systems 
indicates that radio-source structure is independent of group 
environment properties. Our sample covers the full range of typical 
FRI radio morphologies, slightly more than an order of magnitude 
in radio luminosity, and as discussed in the previous section, nearly 
two orders of magnitude in X-ray luminosity. We can therefore 
conclude that, while all of the FRI radio sources in our sample 
have been found to inhabit group-scale environments, the evolution 



of their radio structure appears to be unrelated to the richness of 
the group environment, as parametrized by the X-ray luminosity 
and temperature. The only environmental properties that may be 
relevant to the source evolution are the group pressure and pressure 
distribution. 

5.1.2 Radio morphological class and environment 

To investigate the relationship between radio structure and envi- 
ronment further, we divide the radio sources based on their large- 
scale morphology into the two classes of "plumed" and "bridged" 
sources, as defined by Leahy et al. (1998jf|- By this definition, 
"plumed" FRIs are those in which the majority of lobe material 
is beyond the end of the collimated jet, whereas in "bridged" FRIs, 
the majority of lobe material is between the end of the collimated 
jet and the nucleus. Fig. [9] shows radio maps of a typical member 
of each class. Based on this definition, we classified NCG 1044, 
3C 31, NGC 315, NGC 6251, and 3C449 as "plumed" sources, and 
3C296, 3C76.1, 3C66B and NGC 4261 as "bridged" sources. 

We first tested whether the jet power is the controlling param- 
eter in determining whether a source shows a bridged or plumed 
morphology by comparing the 178-MHz radio luminosities of the 
two subsamples. We found no difference in the radio luminos- 
ity distributions of the two classes, and so we conclude that en- 
vironment must play at least some role in determining the large- 
scale radio structure. We note that for this sample the size distri- 

3 http://www.jb.man.ac.uk/atlas/anatomy.html 
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Table 6. Properties of the radio galaxies and their group environments 



Object 


Ll78 


D 


fax 


V 


^500 


Pe 


Pm 


p 

-* ml 


Re 


Pm 


Ereq 


E rg 


NGC315 


23.38 


1000 


0.14 


65.35 


512 


-14.13 


-13.64 


-14.70 


3.7 


11.4 


B2.86lg;g 


52.31 


3C31 


23.96 


307 


0.16 


63.85 


558 


-13.16 


-12.60 


-13.53 


2.3 


8.5 


52-99-0.5 


51.85 


3C66B 


24.34 


296 


0.35 


64.27 


603 


-12.46 


-12.22 


-13.33 


7.4 


12.9 


CO OQ + 0.15 

OZ.ao_ .2 3 


52.65 


NGC 1044 


23.26 


306 


0.30 


63.82 


473 


-13.25 


-12.98 


-13.82 


3.7 


6.9 




51.44 


3C76.1 


24.40 


126 


0.71 


63.98 


418 


-13.12 


-12.58 


-13.21 


1.2 


4.3 


51 99 +0 14 


52.00 


NGC 4261 


23.67 


81 


0.50 


63.00 


545 


-12.51 


-12.24 


-13.04 


3.4 


6.3 


51 q 2 +"'" 8 


51.36 


3C296 


24.16 


212 


0.46 


64.25 


415 


-13.25 


-12.98 


-13.28 


1.1 


2.0 


< 52.11 


51.87 


NGC 6251 


24.08 


1900 


0.24 


66.08 


577 


-14.63 


-13.83 


-15.70 


12.0 


74 


52.80_°;_* 


52.85 


3C449 


23.81 


533 


0.26 


64.19 


436 


-13.17 


-12.80 


-13.81 


4.4 


10.2 




51.99 



Columns: 1. Radio-source or galaxy name; 2. Logarithm of the 178-MHz radio luminosity (W Hz -1 sr — 1 ); 3. Radio source projected physical size 
(kpc); 4. Axial ratio at widest part of radio source; 5. Logarithm of the sum of the volumes of the two lobes (m 3 ), determined by assuming a spherical 
or cylindrical volume based on the extent of the radio emission; 6. rsoo, the radius enclosing a mean density of 500 times p cr it, determined using the 
relation of Arnaud et al. 1 2005) (kpc); 7. Logarithm of the thermal gas pressure at the outer edge of the radio lobes , averaged for the two lobes except 
in the cases of NGC 3 15 and NGC 6251 (Pa); 8. Logarithm of the thermal gas pressure at the midpoint of the radio lobes, averaged for the two lobes 
except in the cases of NGC 315 and NGC 6251 (Pa); 9. Logarithm of the internal equipartition pressure, averaged for the two lobes except in the cases 
of NGC 315 and NGC 6251 (Pa); 10. Average pressure ratio at the lobe endpoints; 11, Average pressure ratio at the lobe midpoint; 12. Logarithm of the 
energy required to produce the observed temperature offset from the RQ Lx /Tx relation (J) (errors take into account uncertainties in the observed and 
predicted temperatures); 13. Logarithm of the energy available from the current episode of radio-galaxy activity (APmV , J). 




Figure 9. Examples of the two morphological classes of FRI radio galaxies, as described in the text. L: a typical "bridged" source, 3C 296, in which the lobe 
material is primarily between the nucleus and the end of the collimated jet; R: a typical "plumed" source, NGC 1044, in which the lobe material is primarily 
beyond the end of the collimated jet. 



butions of the two classes are significantly different (Fig. II lb . with 
plumed sources being typically larger than bridged sources. Results 
of Kolmogorov-Smirnov (KS) tests comparing the X-ray-measured 
group properties are shown in Table [j_ We find no significant dif- 
ferences between the intrinsic group properties of the two subsam- 
ples, but do find a significant difference between the pressure ratios 
(Pent I Pint) for the two subsamples, as discussed in Section[5__ 

We conclude that the current X-ray data do not conclusively 
demonstrate any close relationship between environmental condi- 
tions and radio-source structure. It appears that both plumed and 
bridged sources can be produced in a range of environments, al- 
though the difference in typical sizes (and axial ratios) of the two 
classes, together with the results discussed in the previous section, 
suggests that the two classes may represent an evolutionary se- 
quence in which the environment plays some role in turning bridge 



structures into plumes. An evolutionary model of this sort would 
require a transition from one morphology to the other. A combi- 
nation of buoyancy effects and collapse due t o the end of over - 
pressuring could produce such a transition (e.g. Hardcastle 1999); 
however, while sources exist in which one side is bridged and the 
other plumed (e.g. 3C 66B), objects in a state of transition between 
bridged and plumed structures do not appear to be common in the 
FR-I population, which suggests that the timescale for the transition 
would need to be short, 

5.2 Source dynamics and particle content 

Fig. [8] shows external pressure profiles for the eight sources ob- 
served with XMM-Newton, determined from the most probable 
pressure at each radius based on the MCMC probability distribu- 
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Figure 10. Correlations between radio structure and group gas properties. Top: Size vs. external pressure (left), pressure ratio (Pext /Pint) (middle), and 
pressure gradient (right); bottom: Axial ratio vs. external pressure (left), pressure ratio (Pext / Pint ) (middle), and pressure gradient (right). 
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Physical size (kpc) 

Figure 11. Histograms of radio-source size for the bridged (hatched) and 
plumed (solid black) sources, showing a statistically significant difference 
in their distributions. 



tion as discussed in the previous section, with internal, equipartition 
pressures overplotted. The plotted uncertainty range takes into ac- 
count uncertainties in both the temperature and the surface bright- 
ness model. Internal equipartition pressures were determined for 
each radio source, and for diffe rent components of the so urce where 
appropriate, using the code of lHardcastle et all l ll998al) . The elec- 
tron energy spectrum was normalized using flux densities measured 
from the 1.4-GHz maps of Table [2] with 178-MHz flux densities 



from the 3CRR catalogue fcaingetaljl983h included where avail- 
able. We assumed *y m in = 10, 7 ma z = 10 6 , an energy injection in- 
dex of 8 = 2 (corresponding to radio spectral index a — 0.5) with 
a break in the spectrum in the region of 7 = 5 x 10 3 — 5 x 10 4 as 
appropriate to fit the radio data points. The equipartition pressures 
assume no contribution from protons. 



5.2.1 Particle content and radio morphology 

In agreem ent with previous studies using earlier X-ray ob serva- 
tories (e.g. lMorganti et alJll988l;IWorr all & Birk inshawlHoOOh and 



using XMM-Newton (e.g. Tcroston et~ 20031) . we found that the 
measured external pressures from the hot gas environments were 
in many cases significantly higher than the minimum internal pres- 
sures for the radio lobes (if it is assumed that the sources are in the 
plane of the sky; relaxation of this assumption does not affect our 
conclusions, as discussed below). Fig.[T2](top) shows histograms of 
the pressure ratios (P ex t /Pint) for the 16 lobes in the sample (only 
one lobe is within the XMM-Newton field-of-view for NGC 315 and 
NGC 625 1) showing that the lobes range from being apparently un- 
derpressured by a factor of ~ 70 to approximate pressure balance. 

This is the first time that such a study has been carried out 
for a sample that includes a wide range of observed FR-I radio 
morphologies, and so a primary aim of our work was to inves- 
tigate whether the amount of "missing" pressure is related to the 
properties of the radio-source and/or its environmental interaction. 
If we consider two extreme cases, 3C 296, which is apparently 
slightly overpressured, and 3C 449, which is apparently underpres- 
sured by a factor of 10, it is clear that 3C 296 is a typical bridged 
source, whereas 3C 449 is plumed. This morphological difference, 
together with the correlations between radio-source structure and 
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Table 7. Correlations between radio-source and group gas properties, de- 
termined using the Spearman rank test. Column 3 is the value of Student's 
t for the given comparison, and column 4 is the null hypothesis probabil- 
ity, i.e. the probability of obtaining the given value of Student's t from an 
uncorrected dataset. 



Parameter 1 


Parameter 2 


t 


NHP 


Lr 


Lx 


0.088241 


0.466 




T x 


0.044102 


0.483 




(3 


1.396017 


0.103 




Tc 


0.683130 


0.258 






0.883277 


0.203 




-* ext j -* int 


0.447 


0.334 




dP/dr 


0.293695 


0.389 


D 


L x 


0.221249 


0.416 




r x 


0.883277 


0.203 




P 


0.221249 


0.416 




r c 


1.272136 


0.12 




Pext 


2.853810 


0.012 




Pext 1 Pint 


2.263 


0.029 




dPIdr 


1.705606 


0.070 


fax 


Lx 


0.634866 


0.273 




Tx 


1.098087 


0.154 




P 


1.098087 


0.154 




r c 


0.988538 


0.178 




Pext 


1.819606 


0.056 




Pext i Pint 


1.396 


0.103 




dPIdr 


1.506237 


0.0914 


Pint 


Lx 


0.088241 


0.466 




Tx 


0.781661 


0.23 




P 


0.176777 


0.432 




r c 


0.935414 


0.190 




Pext 


3.743997 


0.0036 




dPIdr 


1.326473 


0.113 



pressure ratio noted in the previous section, suggests a model in 
which jet/environment interactions affect the particle and/or energy 
content of the radio lobes; for example, the dominant factor in de- 
termining the ratio of radiating to non-radiating particles may be 
the efficiency of the entrainment process. We would expect that the 
entrainment rate for plumed sources is higher, as a larger fraction 
of the jet surface is in direct contact with the external medium, 
whereas a jet that is embedded in a lobe is not in direct contact 
with the group gas (except perhaps very close to the nucleus), and 
so its entrainment rate might be expected to be lower. If the amount 
of entrained material is sufficient to affect significantly the internal 
energy density of the radio lobes (which is likely to require heat- 
ing of the entrained matter), then we might expect to see a larger 
apparent imbalance in the plumed sources. 

The lower panel of Fig. [T2] shows histograms of the pressure 
ratios for the two categories, which demonstrate that there is in- 
deed a signficant difference in the two populations in the expected 
sense: the bridged sources are typically closer to pressure balance 
at equipartition than the plumed sources. Using a KS test, we can 
exclude the hypothesis that the two subsamples are drawn from the 
same parent population at a confidence level of > 97 per cent. As 
shown in Table[8]this does not appear to be solely due to differences 
in the internal pressures for the two subsamples. We note that the 
most extreme pressure ratio in the "bridged" category is the east- 
ern lobe of 3C66B, which has a ratio of 11.5. The classification 
of this source as bridged is not unambiguous, because the eastern 
jet does not enter the lobe for some distance from the nucleus (see, 



Table 8. Comparison of environmental properties for bridged and plumed 
sources. Column 2 gives the KS D value for a comparison between the 
bridged and plumed subsamples, and column 3 gives the null hypothesis 
probability for the subsamples being drawn from the same parent popula- 
tion. For the pressure ratio comparison each radio lobe was treated sepa- 
rately, whereas the other comparisons are for the source as a whole. E rat i 
is the ratio between E req and E rg (see Table|6}. 



Parameter 


D 


NHP 


Lx 


0.50 


0.37 


T x 


0.40 


0.65 


P 


0.50 


0.37 


r c 


0.60 


0.18 


Pext 


0.47 


0.45 


Pext/ Pint 


0.78 


0.005 


Pint 


0.67 


0.10 


dP/dr 


0.50 


0.37 


Ereq 


0.60 


0.26 


E 


0.43 


0.55 




0.60 


0.26 



e.g. lHardcastle et al.l ll996). All of the other bridged lobes have ra- 
tios < 5, and so the distinction between the two populations is 
considerably more significant if the eastern lobe of 3C 66B is ex- 
cluded (a single parent population can be excluded at the 99.8 per 
cent confidence level). 

This comparison therefore supports the hypothesis that en- 
trainment of the jet's surroundings on scales of tens to hundreds 
of kpc is responsible for the apparent pressure imbalance of FR- 
I radio galaxies. Unfortunately there is no obvious observational 
proxy for jet deceleration or entrainment that can be used to test 
this hypothesis furt her with the current data . However, comparisons 
with the models o f lLaing & Bridlelj2002b!, and in prep.) show that 
plumed sources such as 3C31 have an entrainment rate that in- 
creases with distance from the nucleus, whereas bridged sources 
such as 3C 296 may not. Clearly the entrainment characteristics of 
jets that are directly in contact with the external medium will be 
different from those of jets that propagate inside lobe material. The 
strong correlations between pressure ratio and axial ratio discussed 
in the previous section also provide support for this model. 

5.2.2 Projection effects 

The pressure comparisons discussed here use projected distances, 
and hence the true pressure ratios are likely to be somewhat dif- 
ferent. For a radio source at a small angle to the line of sight, the 
source volume will have been underestimated, leading to an over- 
estimate of Pint- At the same time, the external pressure acting 
at the midpoint and outer edge of the lobe (Pext) will be overes- 
timated, because the radio source will extend further out into the 
atmosphere which is decreasing in pressure with radius. Hence the 
extent to which t he pressure ratio is a ffected depends on the par- 
ticular system. In lCroston et al.1 {2003), we argued that for 3C 449 
an angle of < 11° to the line of sight would be required to obtain 
pressure balance, whereas t he two-sided nature of its jets imply 
6 > 75° jFerettietal]|l999t) . For 6 = 75°, the pressure ratio de- 
creases by only 3 per cent from that obtained for the assumption 
that the source is in the plane of the sky (from 12.0 to 11.6 for 
the northern plume). It is clear that projection alone cannot be re- 
ponsible for the large apparent pressure imbalance in this plumed 
source. The source likely to be most affected by projection effects 
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isNGC 625 1 , which has a one-sided jet on large scales. Hones et al.l 
argue that < 40°, based on parsec-scale limits. We find 
that an angle of < 12° is required for the pressure ratio at the mid- 
point of the western plume of NGC 625 1 to be consistent with the 
median value for the bridged sources (P ex t/Pint = 3.3). Such a 
small angle cannot be ruled out, but given the large projected size 
of this source, it is implausible (an angle of < 12° over the entire 
length of the source leads to a physical size > 10 Mpc - see e.g. 
iBirkinshaw & W orrall 1993). An angle of ~ 40°, consistent with 
the pc-scale constraints, would give a more modest reduction in 
pressure ratio, resulting in P cxt / Pint ~ 29. In addition, estimates 
for the line-of-sight angles for other plumed sources (3C 3 1 and 
NGC 315) from kinematic modelling of the rad io data are not ex- 
treme dLaing & Bridld2002al ; Eking et alj2003) . We conclude that 
the difference between the apparent pressure imbalances of bridged 
and plumed sources is not solely the result of projection effects. 

5.2.3 Implications 

If the model we propose here to explain the apparent pressure 
deficit is correct for this sample of FR-Is, then an obvious question 
is whether it is consistent wit h the pressure offsets observed in clus- 
ter centre radio sources (e .g. lMorganti et al ] |l988tlDunn& Fabianl 
|2004 I Birzan et alj l2004h . which typically do not possess well- 
collimated jets on large scales. Given the amorphous structures of 
many cluster centre sources, suggestive of strong environmental 
influence, it seems plausible that there is considerable mixing of 
cluster gas with radio-lobe plasma (and heating, as required in or- 
der for the cavities to be detec ted) before the bu bbles rise in the 
cluster atmosphere. In addition. iDunn e t al. (2006) argue, based on 
a comparison of the leptonic content of the small-scale jet and the 
energetics of the cluster bubbles, that even in bubbles with large ap- 
parent pressure offsets such as those of the Perseus cluster, the jet is 
initially leptonic, so that the pressure offset on large scales is more 
likely to be due to entrainment rather than a relativistic proton pop- 
ulation, consistent with the picture we present here. Hence, while 
the categorization of group-scale sources into bridged and plumed 
morphologies with different energetics may not be directly applica- 
ble to cluster centre sources, these results may nevertheless provide 
some insight into the solution of the pressure balance problem for 
these systems as well. 

An alternative model could be proposed in which the strength 
of the interaction dictates the level of entrainment. For example, 
sources that are overpressured may be less easily disrupted, so that 
entrainment levels are low, while sources close to pressure balance 
are more susceptible to turbulence and higher entrainment rates. It 
is unclear how such a picture would work in detail; however, our 
data do not allow us to rule out such a model. 

5.3 Environmental impact 

As discussed in Section 15.11 it is clear that the group gas sur- 
rounding FR-I radio galaxies is strongly affected by their presence. 
Earlier XMM-Newton studies of the group environments of FR-Is 
found direct evidence for localized heating of the group gas, and in - 
di rect evidence for mor e generalized heating JCroston et al.1120030 . 
In lCroston etall J2005al) we reported an apparent systematic dif- 
ference in the radio properties of galaxy groups with and without 
central radio sources, in the sense that radio-loud groups appear to 
be systematically hotter for a given luminosity. A recent study com- 
paring REFLE2^JORAS£h^rs_with NVSS has found similar re- 
sults ( Magliocchetti & Briiggen 2007), while a Chandra study of a 



subset of the ICroston et al] J2005ah sample by Ijetha et alj fc007h 
has shown that the gas properties in the innermost regions of the 
group do not appear to be significantly affected by the presence of 
a radio source, so that any heating processes must be most signifi- 
cant on larger scales, comparable to the size of the radio lobes. 

In Fig.[T3]we show the X-ray luminosity-temperature relation 
for the current sample of 9 radio galaxy group environments (filled 
circles). The x and + symbol s indica te the radio-quiet and radio- 
loud samples of ICroston et alj d2005al) . respectively, and the solid 
line indicates the best-fitting relation for the radio-quiet groups. 
The points for the XMM-Newton radio-galaxy sample all lie to 
the high-temperature side of the radio-quiet relation, thus confirm- 
ing the earlier result: groups containing extended radio sources 
scatter to higher tem peratures for a given X-ray luminosity. In 
ICroston et a l. (2005a), we argued, based on a lack of difference 
in the X-ray-to-optical luminosity ratios for the two subsamples, 
that this effect must be due to a temperature increase, rather than 
the luminosity decrease expected. The origin of this effect therefore 
remains unclear. It would be expected that any energy injected into 
the environment would alter both its luminosity and temperature, 
with proportionate increases in both if the gas remains in hydro- 
static equilibrium. This is consistent with the position of the FRI 
groups on the Lx-Tx diagram. It has been argued that the slope of 
the Lx-Tx relation is set mainly by the effects of radiative cooling 
(Voigt & Fabian 2004); the offset associated with radio-loud AGN 
activity that we observe here may therefore have the primary effect 
of increasing the scatter in the Lx-Tx relation, rather than altering 
its slope. 

We checked whether it was possible that the presence of 
extended non-thermal emission in the radio-galaxy environments 
could be biasing our temperature estimates by fitting a combined 
mekal plus power-law model to the groups with the largest tem- 
perature excesses, NGC 1044 and NGC 4261, with the temperature 
fixed at that predicted by the radio-quiet Lx /Tx relation. In both 
cases acceptable fits with power-law indices between 1.6 and 2 
were obtained. In the case of NGC 4261 the fit statistic was con- 
siderably worse for this model (x 2 = 1124 for 626 d.o.f.), whereas 
for NGC 1044 the fit statistic was very similar to that of the mekal- 
only model (x 2 = 160 for 129 d.o.f.). The X-ray luminosities of 
the non-thermal components were ~ 6 x 10 41 erg s _1 in both cases, 
comparable to the total luminosity of the group. We cannot, there- 
fore, rule out the presence of extended non-thermal X-ray emission 
in all of the groups in the sample from spectral data alone. To ex- 
plain the entire "temperature excess" in this way would, however, 
require a large fraction of the measured X-ray luminosity to origi- 
nate in extended non-thermal emission. 

Column 14 of Table [6] lists the energy required to raise the 
temperature of all the group gas within the spherical region enclos- 
ing the radio galaxy by AT, the difference between the observed 
gas temperature in that region and the temperature predicted for 
a group of the same luminosity using the radio-quiet Lx /Tx re- 
lation of Crosto net alj d2005al) . To try to understand whether the 
radio source is capable of producing these effects, we estimated the 
total energy available from the current episode of radio activity as 
4PV (this includes the energy stored in the lobes, which cannot 
yet have had an effect on the environment), using the total lobe vol- 
ume and the pressure at the lobe midpoint. Fig.[T3](right) compares 
the required and available energies for each source. As shown in 
Table [8] there is no significant difference in the ratios of required 
to available energy for the bridged and plumed sources. In five of 
the eight sources with temperature excesses the radio source ap- 
pears to be sufficiently powerful to heat the gas by the required 
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Figure 12. Histograms of the P ex t/ Pint, the ratio of external to internal (equipartition) pressure for the full sample. The left panel is for pressures at the 
midpoint of the radio lobe and the right panel for pressures at the outermost of edge of the radio lobe. The sample is subdivided into bridged (hatched) and 
plumed (solid black) lobees. For all except NGC6251 and NGC 315, both lobes are plotted. 



amount, but in three cases (3C31, NGC 4261 and NGC315) the 
available energy is significantly less than that required. It has been 
argued previously (e.g. Nipoti & Binney 2006) that if AGN radio 
oubursts are an intermittent process, then it is most probable that 
we will observe outbursts whose energy injection rate is lower than 
the average for the system. However, if the higher than expected 
temperatures in these systems are the result of radio-source energy 
injection, then they would be expected to disappear on timescales 
comparable to the sound crossing time, so it is not clear that previ- 
ous outbursts from these AGN can be responsible for the effect we 
observe. Another possibility, which is difficult to rule out, is that 
radio-loud AGN prefer to inhabit hotter than average group envi- 
ronments. 



6 CONCLUSIONS 

We have carried out the first XMM-Newton study of a sample of 
low-power (FR-I) radio galaxies and their environments. We find 
the following results: 

• The radio galaxies inhabit group-scale hot-gas environments 
with luminosities ranging from 2 x 10 41 erg s _1 to 10 43 erg s _1 , 
consistent with earlier studies of FR-I environments. 

• The radio properties of the FR-Is, including radio luminosity, 
source size and axial ratio, are uncorrelated with the richness of the 
group environments, as parametrized by either X-ray luminosity or 
temperature. 

• Marginally significant correlations were found between 
source size and external pressure and pressure gradient, and be- 
tween axial ratio and both external pressure and pressure gradi- 
ent, with long, narrow sources associated with lower pressures and 
steeper gradients; however, this could simply be a selection effect, 
as longer sources (which tend to be narrower) are able to probe the 
outer regions of groups, where the pressure gradient is steeper. 

• In agreement with earlier work, we find that FR-I radio galax- 
ies are acted on by external pressures that are significantly higher 



than their internal, equipartition pressures, with apparent pressure 
imbalances ranging from a factor of 70 times underpressured to 
rough pressure balance. 

• We report for the first time that the apparent pressure imbal- 
ance is linked to radio-source morphology: "plumed" FR-I sources 
typically have large pressure deficits, whereas "bridged" FR-Is are 
closer to pressure balance. We interpret this result as evidence that 
"plumed" sources have a higher ratio of non-radiating to radiating 
particles than "bridged" sources, which is likely to be caused by the 
higher entrainment rate expected for "plumed" sources. 

• We find that the temperatures of the FR-I group environments 
are typically significantly higher than predicted for their luminos- 
ity, which supports our earlier claim for a temperature excess in 
radio-loud groups. In 5/8 sources, the current radio source may be 
sufficiently powerful to have produced this excess via heating; how- 
ever, in three cases the current source is too weak. We conclude 
that radio-source heating remains the most plausible explanation 
for this result; however, we cannot rule out an alternative explana- 
tion in which hotter groups are more favourable for the generation 
of an FR-I radio source. 
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Figure 13. Left: The Lx /Tx relat ion for the XMM-New ton radio-galaxy environments sample (filled circles) compared with the radio-quiet (x symbols) and 
radio-loud (+ symbols) samples of Croston et al. ( 2005a). Solid line indicates the best-fitting relation for radio-quiet groups. Right: The energy required to 
produce the observed temperature offsets (CAT) vs. energy available from the current radio source (APV). Solid line is the line of equality. 
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